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Summary. We discuss the origin of physical properties of globular cluster systems 
(GCSs) in galaxies in terms of galaxy formation and evolution processes. Based on 
numerical simulations of dynamical evolution of GCSs in galaxies, we particularly 
discuss (1) the origin of radial density profiles of GCSs, (2) kinematics of GCSs in 
elliptical galaxies, (3) transformation from nucleated dwarf galaxies into GCs (e.g., 
omega Centauri), and (4) the origin of GCSs in the Large Magellanic Cloud (LMC). 



1 Numerical archeology 

Based on penetrative analysis of metal-poor halo stars and globular clusters 
(GCs) in the Galaxy, two canonical Galaxy formation scenarios - the mono- 
lithic collapse scenario and the accretion/merging one j^j- were proposed, 
that have long been influential for later observational and theoretical studies 
of disk and elliptical galaxies. Although observational studies of stellar halos in 
galaxies beyond the Local Group of galaxies have just recently started reveal- 
ing structural and chemical properties of the halos jSJ@], physical properties 
of globular cluster systems (GCSs) in these galaxies have long been investi- 
gated in much more details Wide- field imaging and spectroscopic studies 
with large ground-based telescopes (e.g., Keck 10m) have recently revealed 
GCS structures and kinematics in galaxies with different Hubble types ^6 j. 
Furthermore a growing number of theoretical/numerical studies have recently 
been accumulated which have investigated dynamical and chemical properties 
of GCs and GCSs based on admittedly realistic and self-consistent models of 
GC formation during galaxy formation and evolution 0[H]- In this review 
paper, we therefore try to derive physical meanings from the selected four ob- 
served properties of GCs and GCSs by comparing our numerical simulations 
of GC and GCS formation with latest observations. 
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Fig. 1. Dependences of projected number distributions of stars (thin) and GCs 
(thick) in merger remnants (i.e., elliptical galaxies) on the total number of major 
merger events (./V m ) which an elliptical experienced during its formation (BF06). For 
clarity, the density distributions are normalized to their central values. Thin dotted 
lines represent power-law slopes (a) of a = —2.5, —2.0, —1.5, and —1.0. Note that 
the density profiles of GCSs become flatter for larger N m , i.e. more mergers. 

2 GCS density profiles 

It has long been known that the radial density profiles of GCSs in elliptical 
galaxies (Es) vary with the total luminosities of their host galaxies 9; . If the 
projected GCS density profiles in Es with y-band absolute magnitudes of My 
are fitted to the power-law ones like E gc oc i? Qs ° , where R is the distance from 
the center of the host galaxy of a GCS, the power-law index a gc is smaller 
(i.e., the profiles are steeper) for larger My (i.e., fainter Es). Although two 
physical mechanisms - GC destruction by galactic tidal fields ^01 EH and 
dynamics of galaxy merging (Bekki & Forbes 2006; BF06; ^21) ~ have been 
so far proposed for the origin of GCS density profiles, we focus on the latter 
case in this paper. 

BF06 numerically investigated the structural properties of GCSs in Es 
formed from a sequence of major dissipationless galaxy merging and thereby 
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found that the radial density profiles of GCSs in Es become progressively 
flatter as the galaxies experience more major merger events (See Figure 1). 
The simulated profiles of GCSs are found to be well described as power-laws 
with a gc ranging from —2.0 to —1.0 in Es. They are flatter than, and linearly 
proportional to, the slopes (a s ) of the stellar density profiles. By applying a 
reasonable scaling relation between luminosities and sizes of galaxies to the 
simulation results, BF06 showed that a gc ~ — 0.36Afy — 9.2, r c w — 1.85Afy, 
and a gc « 0.93a s . These correlations between GCS profiles and their host 
galaxy luminosities are consistent reasonably well observations which 
suggests that the origin of structural non-homology of GCSs in Es can be 
understood in terms of the growth of Es via major dissipationless galaxy 
merging. 
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Fig. 2. Correlations between (F m /cro) M p C estimated for R < 2R C and those for 
R < 6R C for metal-poor clusters (MPCs). These correlations are derived from 18 
results of 6 major merger models with three different projections 18 . 



3 GCS kinematics in E/SOs 

Recent observations on GCS kinematics in E/SOs have revealed that GCS 
kinematics can be quite diverse: GCSs in some galaxies like M87 [TJ NGC 
4472 [T3] and NGC 5128 [THj show rotation whereas those in some galaxies 
such as NGC 1399 [T7] do not. Bekki et al. (2005; B05; 18 ) first tried to 
understand the observed diversity in GCS kinematics by numerically investi- 
gating GCS kinematics of E/SOs formed from major/minor galaxy merging. 
B05 demonstrated that both metal-poor cluster (MPCs) and metal-rich ones 
(MRCs) in Es formed from major mergers can exhibit significant rotation at 
large radii (^20 kpc) due to the conversion of initial orbital angular momen- 
tum into intrinsic angular momentum of the remnant. 
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Based on a wide parameter study of galaxy mergers, B05 found that MPCs 
show higher central velocity dispersions than MRCs for most major merger 
models. V m /ao (where V m and a®, are the GCS maximum rotational velocity 
and central velocity dispersion, respectively) ranges from 0.2- 1.0 and 0.1-0.9 
for the MPCs and MRCs respectively, within 6R e for the remnant elliptical. 
Figure 2 shows an interesting result that does not depend on merger param- 
eters: Vm/co of MPCs within 6R e are greater than those of MPCs within 
2R e . 

B05 also revealed the alignment of the major axes in 2D distributions 
between stars, GCs, and dark matter halos in the simulated Es. The aligned 
major axis between stars, GCs, and dark matter appears to be one of the 
principal characteristics of Es formed by major merging, which implies that 
observational studies on 2D distributions of GCSs in Es can tell us about the 
shapes of their host dark matter halos. We also showed in this meeting that the 
total masses of E/SOs estimated from the GCS kinematics can be much closer 
to the real masses than those from the PNe systems owing to less anisotropic 
velocity dispersion in GCSs, in particular, for face-on SOs. This suggests that 
(1) GCSs are better mass-estimators in E/SOs and (2) kinematical data sets of 
PNe systems in E/SOs ^H) should be more carefully interpreted for the total 
masses of E/SOs. Although these simulations are not based on cosmological 
N-body ones, these results may well provide new insight on the origin of the 
observed diversity in GCS kinematics of E/SOs. 
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Fig. 3. The simulated age distributions of field stars (dotted) and clusters (solid) 
in the LMC at the present epoch 1251 . For convenience, the normalized fraction of 
stars in each age bin is shown. 
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4 The age gap problem in the LMC 



Possible candidates of young and metal-rich GCs were discovered in inter- 
acting and merging galaxies [20], and physical properties of these GCs have 
been discussed in different contexts of galaxy and GC formation, such as the 
observed color bimodality in GCSs in Es [Hj , the birth rate of young GCs as a 
function of time in M82 j2J and the age-metallicity relation of GCs in nearby 
interacting galaxies like the LMC and the SMC (e.g., Bekki et al. 2004; B04; 
|22p. In this paper, we focus on the LMC's GCS in which nearly all GCs are 
cither very old (~ 13 Gyr) or younger than 3 — 4 Gyr - the "age-gap" problem 

EHE1 

B04 and Bekki & Chiba (2005; challenged this age gap problem by 
investigating chcmodynamical evolution of the LMC interacting both with the 
Galaxy and the SMC for a long time scale (~ 9 Gyr). They found that the 
first close encounter between the LMC and the SMC about 4 Gyr ago was the 
beginning of a period of strong tidal interaction witch likely induced dramatic 
gas cloud collisions, leading to an enhancement of the GC formation which 
has been sustained by strong tidal interactions to the present day. Figure 3, 
showing the simulated age distributions of field stars and GCs in a model, 
reveals that GC formation can be reactivated about 3 — 4 Gyr ago, when 
the LMC can start its strong tidal interaction with the SMC. These results 
imply that the origin of the age gap can be closely associated with interaction 
histories of the LMC, the SMC, and the Galaxy. 

5 Very massive star clusters 

Very massive star clusters (VMSCs) such as uj Centauri, ultra-compact dwarfs 
(UCDs) , and massive nuclear star clusters have unique characteristics that are 
quite different from those of "normal" GCs. For example, UCDs discovered 
in the Fornax and the Virgo clusters of galaxies j2E]|2Zj|2Hl have intrinsic 
sizes of less than 100 pc, and have absolute B— band magnitudes ranging from 
— 13 to —11 mag, which is more than 2 magnitudes brighter than the most 
massive GC in the Galaxy (i.e., u> Cen). Two physical mechanisms for the 
VMSC formation have been so far proposed: The "Galaxy threshing" scenario 
|2"9~] [3T)| in which VMSCs originate from nuclei of nucleated dwarf galaxies and 
the "cluster merging" one in which VMSCs are formed from merging of smaller 
star clusters in tidal tails of merging galaxies .31 . Although more details on 
physical properties of VMSCs have been recently revealed HUES], ^ remam s 
unclear which of the two scenarios is more convincing for the origin of VMSCs. 

Here we discuss what we can learn from VMSC properties if they were 
previously nuclei of nucleated galaxies. Both dissipationless 34 35 and dis- 
sipative 36 formation scenario of stellar galactic nuclei have provided some 
interesting predictions on nuclear properties of galaxies |37] ■ Figure 4 shows 
the scaling relations between different physical parameters of merger remnants 
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Fig. 4. Correlations of structural and kinematical parameters with My (V— band 
absolute magnitude) for the VMSCs in 40 models [22]. Projected central velocity 
dispersion (ao; upper left), half-light-averaged surface brightness (I e ; upper right), 
effective radius (Re.; lower left), and central surface brightness (iio; lower right) 
are plotted against My. Here the central surface brightness 7io is expressed as 
O.lL/ir / Rio 2 , where L, Rio are the total luminosity of a VMSC and the radius 
within which 10% of L is included, respectively. The best fit scaling relation for the 
VMSCs is derived for each panel using the least square fitting method and described 
as a dotted line with the derived relation (e.g., ao oc L ' 31 ). 



of star clusters initially with the observed GC scaling relations ■ The fact 
that the simulated relations deviate from the GC's ones implies that the scaling 
relations of VMSCs can be used for understanding whether the stellar nuclei 
can be formed from merging of many smaller clusters in the central regions of 
galaxies |39) . The dissipative nucleus formation model |37j has predicted the 
spread of ages and metallicities in stellar populations of stellar galactic nuclei 
and accordingly can be discussed in the context of the observed spread of ages 
and metallicities in lj Cen [6'2\ . Thus dynamical and chemical properties of 
VMSCs can tell us about nucleus formation histories in galaxies. 
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6 Future works: Hierarchical galaxy formation and GCSs 

Thus, structural, kinematical, and chemical properties of GCSs in galaxies 
have fossil information on dynamics of major/minor galaxy merging (e.g, an- 
gular momentum redistribution processes in merging), interaction histories of 
galaxies, and the formation histories of stellar galactic nuclei in dwarfs. Pre- 
vious theoretical/numerical studies however did not discuss so extensively the 
observed correlations between GCS properties and their host ones in the 
context of a hierarchical clustering scenario of galaxy formation. Several au- 
thors just recently have started their investigation on the GCS-host relations 
based on semi-analytic models |40| and high-resolution numerical simulations 
in ylCDM models A number of observed GCS-host relations, such 

as the positive correlation between GCS metallicities and their host galaxy 
luminosities [5] have not been clearly explained by any galaxy formation sce- 
narios. Since these GCS-host relations may have profound physical meanings 
on galaxy formation and evolution, it is doubtlessly worthwhile for future nu- 
merical simulations of GCS formation to explore the origin of these relations. 
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